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SECTION  1 
INTRODUCTION 


When  trying  to  simulate  SGEMP  in  a vacuum  tank,  cavity  resonances 
are  almost  always  excited.  Exciting  these  resonances  may  or  maynot  affect 
measurements  made  on  objects  within  the  tank.  In  this  report  we  will  be 
concerned  with  calculating  the  amplitudes  of  cylindrical  cavity  resonances 
under  a special  set  of  circumstances. 

We  shall  consider  circumstances  which  approximate  severe  space- 
charge  limiting.  The  objective  of  these  considerations  will  be  to  obtain 
an  understanding  of  the  excitation  amplitudes,  in  general,  as  functions  of 
the  time  and  length  parameters  involved  and  also  to  provide  some  specific, 
relevant,  theoretical  computation  for  the  experiments  performed  by  Mission 
Research  Corporation  at  Physics  International. 

In  a simulation  of  SGBIP,  in  a tank,  photons  are  allowed  to  impinge 
on  an  object  within  the  tank.  The  photons  cause  photoelectrons  to  be  ejected 
from  the  object  and  the  tank  walls.  These  photoelectrons  form  sources  of 
electromagnetic  fields;  the  tank  walls  and  object  are  the  boundaries  for 
the  electromagnetic  fields.  In  reality  the  tank  and  the  object  are  the 
system  which  has  resonant  modes. 

In  our  calculations  we  take  the  point  of  view  that  the  source  cur- 
rents are  axially  symmetric  and  that  the  presence  of  an  object  does  not 
greatly  effect  the  modes  of  the  cavity.®  We  also  assume  that  the  sources 
can  be  separated  into  the  product  of  a function  depending  only  upon  spatial 
variables  and  a function  depending  only  upon  time.  In  Section  2.1  we  obtain 
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a general  formula  for  the  modal  amplitude  of  the  magnetic  field  under  these 
assumptions. 

Under  conditions  of  severe  space  charge  limiting  the  sources 
reduce  to  electric  dipole  moments  in  space.  If  the  test  object  is  a disk 
on  which  the  dipole  moment  is  uniform,  the  source  of  fields  exists  only  at 
the  edge  of  the  disk  where  the  dipole  moment  is  discontinuous.  In  Section 
2.2,  using  the  formula  developed  in  Section  2.1,  we  derive  a formula  for 
the  modal  amplitude,  excited  by  one  point  of  discortinuity,  after  the  source 
currents  have  ceased.  Modal  excitation  due  to  an  arbitrarily  shaped  object 
could, however, be  approximated  as  a sum  of  point  discontinuities. 

A discussion  of  the  consequences  stemming  from  the  general  form 
of  the  equations  derived  in  Section  2.2,  together  with  numerical  evaluations 
is  presented  in  Section  3.  The  tank  dimensions  used  in  these  numerical 
evaluations  are:  tank  radius  equal  to  2 meters  and  tank  length  equal  to  6 
meters.  These  dimensions  are  loughly  equivalent  to  the  Physics  International 
simulation  tank  used  with  Owl  11.  Table  1 presents  some  of  the  periods 
for  the  modes  of  this  tank.  Graphs  which  should  simplify  modal  excitation 
<3StimaLes,  for  an  arbitrary  tank,  are  also  presented  in  Section  3. 

The  results  and  conclusions  are  summarized  in  Section  4. 
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SECTION  2 
THEORY 


2.1  GENERAL  RELATION 


In  this  section  we  derive  the  general  relation  which  expresses  the 
amplitude  of  a cylindrical  cavity  mode  in  terms  of  the  position,  magnitude 
and  time  history  of  the  driving  currents.  We  begin  with  Maxwell's  equations 
for  an  axisymmetric  situation.  By  a simple  manipulation  of  these  equations, 
the  equations  describing  the  fields,  in  cylindrical  coordinates,  become: 


where 


c^  3t^ 


4tt  / 3 3 T \ 

® ^ c 1 3r  '^z  ~ 3z  rj  ’ 


(2-1) 


(2-2) 


, 3E 
^ r 

c 3t 


^■n  3B 

c '^r  ” 3z 


(2-3) 


and  also 


1 

J z_ 

c 3t 


ii  I + 1 frB'l 


(2-4) 


In  the  above  equation  - subscript  z or  r refers  to  the  z or  r cylindrical 
coordinates  respectively;  E and  B refer  to  the  electric  and  magnetic  fields 
respectively;  and  J refers  to  the  spatial  current  density.  If  Equations  2-1 
through  2-4  are  to  be  satisfied  in  a conducting  cavity  (E  fields  parallel 
to  the  conducting  surfaces  are  zero),  O^rsR,  0<zsL  then 


and 


Hr.z.t)  . f t 

n"l  p*l  ^ ' ' 


:2-5) 


(2-6) 


where  are  Bessel  functions  of  order  one  and 

R L 


= 4(R2uJ(X„))-‘/rdr/d.Wi(-ir)co.(Ep.)  . (2- 


2-7) 


0 0 
R L 


g^p(t)  . 4(R2uJ(X^))’^/rdr/d2gJj(-^r)cos(2p-)  . (2-8) 


0 0 


_3_ 

3r 


(2-9) 


r»R 


In  Equations  2-7  and  2-8  we  have  used  the  normalisation  relationship  for  the 
Bessel  functions  I-IO,  where  KR  in  relation  I- JO  is  replaced  by  the  root 
The  roots,  are  found  from  the  boundary  condition  Equation  2-9. 


We  now  need  to  find  Bpp(t)  in  terms  of  gj^p(t)  to  solve  the  problem. 
Substituting  Equations  2-5  and  2-6  into  Equation  2-1  we  have 


3^B 


np  . 2 

■s*-  10 


3t 


2 np“np 


= - c' 


^p 


(2-10) 


where 


2 2 
0)  = c 

np 


[(t)  *(^)]- 

Equp.tion  2-10  can  be  solved  to  yield 


(2-11) 
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t 


np 


(t)  = Jc  (t*)  sin  0)  (f-t)df 


(2-12) 


np 


where  we  assume  that  g(t)  ■=  0 for  t < 0.  Combining  Equations  2-2,  2-8  and 
2-12  we  obtain  the  equation  we  are  seeking: 

R LJj(X^)u)^p  0 0 0 ‘■ 

•^r)]  * (2-13) 


2.2  BOUNDARY  LAYER  AND  TRIANGULAR  TIME 
HISTORY  PARTICULARIIATION 


Having  obtained  the  general  relation  describing  the  amplitude  of 
a cylindrical  cavity  mode  as  a function  of  time,  for  a general  driving 
current,  we  calculate  the  amplitude  for  a circumstance  which  represents  the 
setting  up  of  a boundary  layer  on  a disk.  The  disk  has  a radius  r^  and  is 
located  at  the  position  2^.  The  axis  of  the  disk  corresponds  to  the  axia  of 
the  cavity.  The  spatial  currents  are  idealized  as 

= 0 , (2-14) 


6(2-ZQ)H(rQ-r)6(t)  , 


(2-15) 


where  we  first  use  a very  simple  time  history, 

S(t)  = -7t  0<t<T, 

T 

S{t)  =0  0 > t > T , 


(2-16) 


In  Equation  2-15  X represents  a length  characteristic  of  the  boundary  layer 
thickness;  represents  a characteristic  surface  charge  density  for  the 
boundary  layer  and  T represents  the  time  it  takes  to  set  up  the  boundary 
layer;  6 and  H are  a delta  function  and  a step  function  respectively,  he 
next  substitute  Equations  2-14  through  2-16  into  Equation  2-13  and  perform 
the  specified  integrations. 
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The  r and  z integrations  are  as  follows: 
L R 


/dz/dr  r(^  J^|cos(pTTz/L)Jj  = - aQXrQS(t)Jj(^  rQ]cos(pTTZQ/L) 


0 0 


(2-17) 


The  integration  over  time  is,  for  t > T 

T T 


/S(f)sin((ij^p(t'-t))df  = - ^ cos(w^p(f-t)dt' 


0 


•-•‘p  0 

T 


np 


np  0 


i-  1*-^-  S(t')sin(w  (t’-t))df  . 


0)  0 3t 

np  ^ 


y (sin((jjj^pt)-sin(uJj^p(t-T))  , (2-18) 


since 


-^S(t)  = ~ (5(0)  - 6(t-T))  . (2-19) 

3t  T 


We  further  reduce  Equation  2-18  by  using  the  formula  for  the  addition  of 
sine  functions 

sin(Wj^pt)-sin(a)^p(t-T))  = 2cos(Wj^p(t-T/2))sin((;^^pT/2)  . (2-20) 

Substituting  Equation  2-20  into  2-18  we  have 

T 2 

/S(t')sin(w  (t'-t))dt'  = (jif-)  COSO)  (t-T/2)(sin(w^pT/2))  . (2-21) 

0 "^P 


Using  Equations  2-17  and  2-21  in  2-13  we  have 


B (t) 
np^ 


64ircp»Ar 
1 ^ n-^  np 


Using 


01  = 27T/t 

np  np  • 


(2-23) 


where  x^p  is  the  period  o£  the  mode  described  by  th“  subscripts  np,  we  can 
express  Equation  2-22  as 


/VVo^\ 


(cos(p  Zq/L))cos(w  (t-T/2))  (2-24) 


If  we  had  defined  S(t)  by  a pulse  which  rises  in  T and  decays  to 


zero  in  2T, 


t 0 < t < T , 


S(t)  = -y  (2T-t)  T < t < 2T  , 
T 


S(t)  = 0 2T  < t , 


(2-25) 


we  would  have  obtained  Equation  2-24  with  sin  (irT/x^p)  replaced  by  sin'^ 
(TiT/x^p)  and  cos(w^p(t-T/2))  replaced  by -sin (w^^p (T-t) ) . That  is,  with  this 
new  time  history. 


!).  ijiji 


np 


= “o(;^)a)(-^)(f)(¥)'-"^(^) 

/JlCX^ro/R)\ 

( ,2,,,  )(cos(pTTZQ./L))(sinw  (t-TO) 


(2-26} 


The  pulse  defined  by  Equation  2-25  approximates  the  curve  describing  the  time 
derivative  of  the  dipole  moment  of  a linear  times  exponential  energy  distri- 
bution with  a linearly  rising  pulse. ^ [We  will  call  the  pulse  described  by 
Equations  2-16,  pulse  1 and  that  described  by  Equation  2-25, pulse  2.]  Since 
pulse  2 approximates  a real  physical  situation.  Equation  2-26  will  be  referred 
to  more  often  than  Equation  2-24.  Pulse  1 is  designed  to  indicate  the  ef- 
fects of  a very  sharp,  more  pathological  pulse.  The  actual  pulse  in  Reference 
2,  page  29,  Figure  23,  does  not  really  end  as  rapidly  as  the  pulse  described 
by  Equations  2-25.  Electrons  escaping  the  boundary  layer  tend  to  lengthen 
the  time  over  which  the  dipole  moment  changes.  In  our  calculations  we  are 
not  considering  these  electrons, since  their  effect  may  not  be  adequately 
treated  by  considering  the  boundary  layer  thickness  to  be  small  compared  to 
the  dimensions  of  the  emitting  object.  In  addition,  real  photon  pulses  may 
linearly  rise  for  only  a few  boundary  layer  rise  times;  the  photon  pulse  in 
Reference  2 rises  continually. 


The  required  Equation  2-24  expresses  the  amplitude  of  transverse 
magnetic  modes  described  by  the  indices  n and  p in  terms  of  the  sources; 
OpASCt)  is  the  rate  of  change,  with  respect  to  time,  of  the  dipole  moment 
per  unit  area.  Equations  2-18  and  2-19  show  that  the  amplitude  of  the  mode 
for  these  S(t),  can  be  directly  expressible  in  terms  of  the  third  derivative, 
with  respect  to  rime,  of  the  dipole  moment  per  unit  area.  In  other  words 
a contribution  is  made  to  the  amplitude  of  a mode  every  time  the  slope  of 
the  S(t)  curve  changes.  These  changes  can  cancel  each  other.  Such  is  the 


case  when  T 


(see  Equation  2-24  or  2-26) 
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SECTION  3 


-1 


RESULTS 


3.1  SOME  GENERAL  CONSIDERATIONS 

Equation  2-24  and  2-26  describe  the  amplitude  of  the  n,p  mode  after 
the  driving  current  pulse  is  over.  These  expressions  involve  a factor  which 
depends  upon  the  roots  of  the  Bessel  function  Jp  and  a factor  which 
depends  upon  the  ratio  of  the  rise  time  to  the  period  of  the  mode.  These 
factors  are  plotted  in  Figures  1 and  2 respectively.  Figures  1 and  2 can 
be  used  to  estimate  cavity  modal  r'?sponses.  Figure  2 contains  plots  for 
the  two  S(t)  discussed  in  Section  2.  For  small  values  of  T/r  the  curve 
describing  pulse  1 goes  as  t'’ (Tjjp/T)  and  the  curve  describing  pulse  2 goes 
as  'iT^  (Xj^p/T) . 

In  Section  3.2  we  will  discuss  the  distribution  of  amplitudes  of  the 
n,p  modes  with  respect  to  the  amplitude  of  the  lowest  mode,  n = l,  p = 0.  We  will 
now  discuss  how  the  amplitude  of  the  modes  depend  upon  the  various  parameters 
in  Equation  2-24  and  2-26.  It  is  clear  from  both  these  equations  that 
smaller  values  of  T in  general  give  rise  to  larger  amplitudes,  however  in 
the  limit  that  ttT  « the  amplitude  for  the  np  mode  with  pulse  2 is 
independent  of  T.  For  the  sharper  pulse,  pulse  1,  the  amplitude  increases 
as  T"^  for  T ->■  0.  From  Figure  2,  for  n = 1,  the  amplitude  of  the  lowest 
mode  is  smallest  for  smaller  values  of  r-/R.  In  the  limit  that  X r^/R  0 

2 V n U 

the  amplitude  of  np  mode  goes  as  (r^/R)  (see  Equation  2-24  and  2-26) . 

Keeping  the  points  of  discontinuity  close  to  the  axis  for  a given  tank  will 
therefoie  reduce  the  amplitude  of  the  lower  n modes  or  alternatively 

11 


increasing  the  radius  of  a tank,  R,  over  the  object  size,  r^,  will  reduce 
the  anplitudes  of  the  inodes  like  (r^/R)^. 


We  will  now  consider  how  the  amplitudes  depend  upon  the  physics 
of  the  boundary  layer.  If  we  imagine  that  the  modes  are  stimulated  by  a 
linearly  rising  photon  pulse  with  a linear  times  exponential  emission, 
electron  energy  distribution, then  from  Reference  2 


OqX 


- statcoul 
.417  £ 


cm 


(3-1) 


and 


gl/2  1/3  2/3 

T . 1.28 


sec 


(3-2) 


where  ¥ is  the  e-folding  energy  of  the  distribution  (in  keV) , Y is  the 
material  yield  (elact/cal),  f is  fraction  of  the  energy  in  the  rise  of  the 
photon  ^ulse,  tj^  is  the  rise  time  of  the  photon  pulse  (in  sec).  The 

amplitude  of  the  cavity  modes  depends  upon  the  factor  o^X/T  (times  a 
trigonometric  function)  which  from  Equations  3-1  and  3-2  is 


°o'' 


(t 


I 


(3-3) 


The  amplitude  of  the  modal  response  is  therefore  mere  sensitive  to  the  rise 

time  of  the  photon  pulse  than  it  is  to  the  other  factors  for  T s,  T^^p.  For 

T '»  T the  amplitude  of  the  modal  response  is  proportional  E , with  pulse 

_ 5/6  2 1/3 

2;  with  pulse  1 the  amplitude  is  proportional  to  E (Y4)f/tj^)  . In  other 

words  for  very  short  boundary  layer  rise  times  the  average  energy  of  the 
electrons  is  the  more  important  quantity  for  stimulating  the  lower  frequencies. 
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3.2  FIELD  DISTRIBUTION  AMONG  MODES 


It  is  important  to  understand  what  the  relative  amplitude  of  the 
modes  are  under  various  conditions,  especially  if  one  is  interested  in  damp- 
ing these  modes.  Tables  2-6  show  the  variation  of  modal  distribution, 
normalized  to  the  lowest  mode,  as  T is  varied  from  1 NS  to  9 NS,  for  pulse  2. 
Mode  amplitudes  are  given  for  the  modes  1 S n £ 5 and  0 £ p <9.  The  tank 
parameters  are:  R •=  2 m and  L = 6 m.  The  position  of  the  discontinuity  is 
at  Tq  e .5  m,  Zq  « 1.5  m.  These  spatial  parameters  correspond  to  those  of 
the  "Disk"  experiment  performed  by  Mission  Research  Corporation  (MRC)  at 
Physics  International  in  August  of  1976.  Table  1 gives  the  periods  of  the 
n,p  modes  for  the  given  tank  parameters. 

It  is  worthwhile  noting  from  Tables  2-6,  that  for  a boundary  layer 
rise  time  which  is  smaller  than  2.5  NS  some  modes  have  amplitudes  greater 
than  the  first  mode.  Combining  the  amplitudes  of  Tables  2-6  with  the  phase 
factors  given  in  Equation  2-26  and  inserting  the  result  in  Equation  2-5 
we  obtain  Figures  3-7  and  8-12.  Figures  3-7  are  plots  of  the  normalized 
magnetic  field  at  r = .35  m,  z = 1.5  m and  Figures  8-12  are  plots  of  the 
normalized  B at  r = 2 m and  z = 3.3  m.  These  positions  correspond,  roughly, 
to  two  sensor  positions  in  the  MRC  "Disk"  experiment,  one  on  the  experimental 
object  and  the  other  on  the  tank  wall.  As  the  rise  time  of  the  boundary 
layer  increases  the  lowest  mode  becomes  more  obvious  as  part  of  the  "noise." 
The  lowest  mode  appears  to  be  more  obvious  at  the  tank  wall  position  than  at 
the  position  of  the  object. 

We  now  examine  the  distribution  of  modes  from  a more  general  point 
of  view.  If  represents  the  amplitude  of  the  np  mode, then  from  Equation 
2-26  for  pulse  2, we  have 
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Figure  10.  Sum  of  modes  vs.  time,  T = 5.0  NS 


Figure  12.  Sum  of  modes  vs.  time 
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Since  > T it  is  clear  from  Equation  3-4  that  smaller  values  of  T 
10  np 

favor  larger  relative  amplitudes  for  the  n,p  modes.  Smaller  values  of  i'q/R* 

in  general  for  a given  n,  also  favor  larger  relative  amplitudes  for  the  n,p 

modes.  For  r^/R  = 1,  the  factor  dependent  upon  r^/R  in  Equation  3-4  is  1. 

For  Cq/R  0,  the  factor  dependent  upon  r^/R  in  Equation  3-4  is  equal  to 

Xn/Xi*  where  X^  > X^^.  (The  factor  dependent  upon  r^/R  goes  through  zero 

as  Tq/R  is  varied  between  1 and  zero  but  the  trend  is  for  the  extrema  of 

the  factor  to  increase  as  0*  I't  should  also  be  noted  that 

2 2 

J^(xp/Jj(X^)  is  an  increasing  function  of  n so  that  for  T « and 
\ Tq/R  « 1 Equation  3-4  becomes,  with  the  help  of  Equation  2-11 


Ji(xp 

J?(X  ) 
I'-  n^ 


(3-5) 


For  p 


0 or  R/L  <<  1 we  have 


B 

^ s 

^10 


J?(X  ) 

1 ^ n"^ 


(3-6) 


that  is  under  extreme  conditions  the  higher  n modes  actually  dominate,  in 
a way  specified  by  Equation  3-6,  over  the  lowest  mode,  since  'ij^(Xj^)  ^ Jj^(Xp. 
The  relative  amplitude  dependence  on  z position  is  given  simply  by  the 
cosine  function  in  Equation  3-4. 

For  a particular  position  within  the  tank  the  mode  amplitudes 
continue  to  give  rise  to  a total  magnetic  field  through  Equation  2-5.  The 
extrema  of  Jj(y)  are  a decreasing  function  of  y so  the  general  trend  is 
for  modes  with  larger  n's  to  contribute  less  to  the  field  as  the  tank  wall 
is  approached. 
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3.3  APPROXIMATE  MODAL  MAGNITUDES 


We  calculate  the  amplitude  of  the  lowest  mode  (n  = 1,  p = 0)  for 
T = 2.5  NS  and  T = 5 NS  where  r^/R  = z^/L  = .25.  If  E = .6  keV  then  from 

Equation  3-1  a^X  = .25  statcoul/cm.  For  T <=  2.5  NS,  from  Table  1 T/Xj^p  ■ 

.143;  for  T = 5.0,  from  Table  T/x^^p  = .286.  Looking  at  Figure  2,  for  the 
pulse  2 curve,  at  these  ratios  we  find  that  the  time  factor  is  about  60  and 
25  respectively.  We  next  look  at  Figure  1 for  the  Bessel  function  factor; 
it  is  about  1.  Inserting  the  required  number  into  Equation  2-24  we  find 
that 

s 2 X 10“^  gauss  , 

for  both  values  of  T.  Looking  at  Figures  3 and  4 we  might  expect  the  maxi- 

-3  -3 

mum  B field  of  the  noise  to  be  about  4 x lO  to  1 x lO  gauss  for  T *=  2.5 

and  T = 5.0  at  the  body  sensor  position.  If  the  rise  time  of  the  boundary 

layer  were  only  one  nanosecond  one  might  expect  a peak  noise  magnetic  field  of 
about  10  X 10  gauss. 
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SECTION  4 


SUMMARY  AND  CONCLUSIONS 


In  Section  2 an  equation  was  derived  which  describes  the  modal 
excitation  of  a cylindrical  cavity  for  a highly  space-charge -limited  situa- 
tion] the  space-charge  dipole  moment  was  uniform  except  at  one  point.  A 
discussion  of  this  equation  was  given  in  Section  3 together  with  some  numeri- 
cal evaluations.  One  of  the  time  histories  discussed  in  Section  3 approxi- 
mated that  of  a detailed  one-dimen.sional  calculation.  The  one-dimensional 
calculation  assumed  a linear  times  exponential  electron  energy  distribution 
and  a linearly  rising  photon  time  history. 

When  the  amplitude  of  the  excited  modes  are  expressed  in  terms  of 
the  boundary  layer  parameters  it  appears  that,  for  boundary  layer  rise  times 
which  are  larger  than  or  roughly  equal  to  a modal  period,  the  rise  time  of 
the  photon  pulse  is  the  most  important  parameter  in  exciting  that  mode.  If 
the  rise  time  of  the  boundary  layer  is  much  smaller  than  a modal  period  the 
average  energy  of  the  ejected  electrons  become  the  predominant  factor  in 
exciting  that  mode. 

Qualitatively  speaking,  shorter  boundary  layer  rise  times  favor 
larger  amplitudes  and  increase  the  coupling  to  higher  modes.  As  the 
rise  time  approaches  zero,  the  distribution  of  amplitudes,  relative  to  the 
lowest  mode,  approaches  specified  values  independent  of  the  rise  time  but 
dependent  upon  the  position  of  the  point  of  discontinuity  of  the  dipole 
layer.  These  conclusions  are  what  one  generally  expects  but  the  equations 
in  the  text  transf jrm  the  general  assertions  into  exact  magnitudes. 
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As  the  rise  time  of  the  boundary  layer  decreases, the  amplitudes 
of  the  modes  become  approximately  independent  of  rise  time  for  one  case 
considered:  that  of  a boundary  layer  whose  electrons  have  a linear  times 
exponential  distribution.  Again,  qualitatively  speaking,  a point  of  dipole 
discontinuity  close  to  the  axis  favors  higher  frequencies  and  smaller 
amplitudes.  The  amplitude  of  the  lower  modes,  for  a point  of  discontinuity 
close  to  the  axis  (or  a small  body  in  a large  tank)  goes  roughly  as  the 
volume  of  the  space-charge  layer  over  the  volume  of  the  tank. 

Determining  the  amplitude  and  distribution  of  modes  is  an  important 
consideration  in  the  simulation  quality  of  any  tank,  especially  if  the  modes 
must  be  damped.  (Higher  frequencies  appear  to  be  more  easily  damped  than 
lower  frequencies.)  Numerical  evaluation  for  spatial  parameters  relevant  to 
the  P.I.  tank,  and  a range  of  boundary  layer  rise  times  indicate  that  the  modal 
excitation  is  not  always  dominated  by  the  lowest  mode.  A damper  grid  for  this 
tank,  to  be  used  under  severe  space-charge  limited  circumstances,  should  be 
designed  to  be  effective  for  a range  of  frequencies  higher  than  the  lowest  modOv 

In  a simulation  tank  electrons  emitted  at  the  walls  of  the  tank  also 
contribute  to  modal  stimulation.  If  these  wall  electrons  are  space  charge 
limited  then  the  equations  in  the  text  of  this  report  can  be  used  to  estimate 
the  amplitude  of  excitation  of  axisymmetric  modes.  The  modal  amplitude  caused 
by  a source  discontinuity  at  the  curved  surface  of  the  tank  is  less  than  for 
the  optimum  point  of  stimulation  for  that  mode  (the  optimum  point  of 
di.scontinuity  would  occur  at  some  radial  point  smaller  than  the  radius  of  the 
tank);  the  distribution  of  modes  for  amplitudes  stimulated  at  the  curved  sur- 
face of  the  tank  tend  to  be  towards  the  lower  frequencies.  The  amplitude  for 
a source  continuity  at  the  flat  surfee  of  the  tank  is  larger  or  equal  to  the 
amplitude  stimulated  by  a discontinuity  elsewhere  in  the  axial  direction. 
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APPENDIX  I 


In  this  appendix  we  derive  the  normalization,  for  Bessel  functions, 
used  in  the  body  of  this  report.  Formulas  for  normalizing  Bessel  functions 
appear  in  ne  literature.  However,  these  formulas  are  not  always  consistant 
from  author  to  author.  This  short  proof  may  satisfy  the  skeptic.  The  proof 
is  designed  to  examine  the  boundary  conditions  which  give  rise  to  the  normal- 
ization formula.  We  begin  by  obtaining  a useful  relation  from  Bessel's  equation. 
Bessel's  equation  is 


r 


(I-l) 


Multiplying  Equation  I-l  by  r and  integrating,  by  parts,  in  the  range 

a < r £ b 


r,  2,  9 T 1 r r 2.2  2,  9 , .2. 

I drr  J J ■=  — (n  J - r J ) ) 
J n 9r  n _,,2  I ^ n ^9r  n'^ 

2K  '-a 


(1-2) 


where  we  have  used  the  fact  that 


2 99  1 o 

/ drJ  = i r 

J n 9r  2 (_  n 


« I-  jn 


(1-3) 


(1-4) 


With  Equation  1-2  we  are  in  a position  to  find  the  normalization  constant. 
It  is  clear  that 
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1 


8 . T ^ T ^ T 

(rJJ  = * r ^ . 


dr  ^ n" 


n 


(1-5) 


t-lultiplying  1-5  by  rJ^  and  integrating  by  parts  over  the  range  a < r < b we 


find  that 


r rJ^dr  = ^ r r^ - /*  drr^ J J , 

J n 2 L n J n 3r  n ' 


(1-6) 


substituting  Equation  1-2  into  1-6  we  find  that 

b 


(1-7) 


For  the  normalization  considered  in  this  report: 
a = 0 , 
b «=  R , 
n «=  1 ; 

and  also  for  the  boundary  condition: 


(1-8) 


3r  1 


Jl(KR) 


(1-9) 


r=R 


we  find  upon  substituting  Equation  1-8  and  then  1-9  into  Equation  1-7  that 
R 

/rJ^dr  = ~ (RJ^(KR))^  . (I-IO) 


0 


Equation  I-IO  is  not  the  same  as  that  found  on  page  73  of  Reference  1,  for 
example. 


! 1 
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AT'iN:  William  L.  Chadsey 

Spire  Corporation 

ATTN:  Roger  G.  Little 

SRI  International 

ATTN:  Library 

Systems,  Science  and  Software,  Inc. 

ATTN:  Andrew  R.  Wilson 
ATTN:  Technical  Library 

TRW  Defense  & Space  Sys.  Group 

ATTN:  Tech.  Info.  Center/S-1930 
2 cy  ATTN:  Robert  M.  Webb,  Rl-2410 
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